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Abstract. The experimental results obtained for the magneto-transport in the 
InGaAs/lnAlAs double quantum wells (DQW) structures of two different shapes of 
wells are reported. The beating-effect occurred in the Shubnikov-de Haas (SdH) 
oscillations was observed for both types of the structures at low temperatures in the 
parallel transport when magnetic field was perpendicular to the layers. An approach 
to the calculation of the Landau levels energies for DQW structures was developed 
and then applied to the analysis and interpretation of the experimental data related to 
the beating-effect. We also argue that in order to account for the observed magneto- 
transport phenomena (SdH and Integer Quantum Hall effect), one should introduce 
two different quasi-Fermi levels characterizing two electron sub-systems regarding 
symmetry properties of their states, symmetric and anti-symmetric ones which are 
not mixed by electron-electron interaction. 
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1. Introduction 

In recent years a great number of publications appeared, where the experimental studies 
of electron transport in Double Quantum Well (DQW) structures were reported. The 
reason for this is the enduring progress in the structure fabrication technology which 
allows to modify the shape and width of QW or barriers, as well as the progress 
in measurement techniques enables to measure the separate conductivity and Hall 
response of each layer. These experiments have demonstrated a number of new and 
interesting phenomena occurring in such structures which however, cannot be observed 
in a structure with a single two-dimensional-gas (2DEG) layer. Generally, all these 
phenomena can be divided into two groups depending on the intensity of interlayer 
interaction. To the first group the phenomena related to the Coulomb drag effect belong, 
while to the second one, the phenomena related to the tunneling effect. Correspondently, 
the phenomena of the first group are accounted for the weak interaction between the 
layers, while the phenomena of the second group are attributed to the strong interaction. 

In order to reveal more clearly drag effect caused by Coulomb interaction, the 
tunneling between two QWs has to be negligible and in order to fulfill this requirement, 
the barrier between two QWs should be adequately wide. M.C. Bronsager et. al[T] 
reported theoretical results concerning the magnetic-field and temperature dependences 
of the trans-resistivity in DQW-layers. Experiments of Ref. 2 confirmed the results 
obtained in Ref. 1 and have shown also the existence of double-peak structure for 
each magneto-resistance peak in the samples characterized by the matched density of 
states in the layers. Another type of experiments were carried out on the samples with 
unmatched density of states[3]. In this work authors demonstrated the existence of 
negative drag in the IQHE-regime; its clear explanation was proposed by J. G. S. Lok 
et. al[3]. The authors of Ref. 4 argued that the electron spin plays decisive role in the 
drag effect. At last, the fractional drag between parallel two-dimensional systems has 
been observed and measured in a regime of strong interlayer correlations [5], where the 
Hall drag resistance was observed to be quantized at h/e^ . 

Another situation takes place if one has two non-identical QWs. In this case some of 
the integer Hall states appear and some others disappear, which depends on alignment 
of electron states in two QWs[6]. 

A lot of new phenomena caused by the tunneling between two 2DEG-layers 
separated by the barrier were also observed in recent years. One of them is the absence 
of some integer quantum Hall states caused by the splitting of electron states into 
symmetric and anti-symmetric ones (so colled SAS-gap, or A^as) in DQW systems [71 [8]. 

In a structure where the barrier is thin enough, the SAS-gap plays decisive role. 
The magnitude of this splitting can be determined directly by means of current- voltage 
characteristics of such structure; it was proven also that the SAS-gap is proportional to 
the magnetic field[9]. This effect was attributed by D. Huang and M.O. Manasreh|10] 
to the screening of electron-electron interaction in the DQW-structure. 

In a last few years many new phenomena were observed in the structures 
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with strongly correlated 2DEG at very low temperatures, among them the Bose 
condensation [1 1 j . Wigner crystallization [121 [IS]) and the long-time nuclear-lattice 
relaxation in DQW[lll [15]. Some authors[l6l [IT] have considered the Shubnikov-de 
Haas (SdH) oscillations in a tilted magnetic field with the strong perpendicular 
and small parallel B\^\ components of the field. The authors of Ref. 17 stated that the 
small B\\- component is responsible for the beating effect in SdH oscillations occurred in 
DQW-structure. If there is no i?x-component and magnetic field is parallel to the layers 
composing DQW, the most important peculiarities of magneto-transport are accounted 
for the crossing of Fermi- surfaces belonging to two QWs[T8l [191 120] and this case is 
completely different from two previous ones. 

Therefore, as for the electrons in strongly correlated DQW-structures is concerned, 
it is possible to treat them as two separate sub-systems, one of them as belonging to 
symmetric states, while another one as to belonging to the anti-symmetric ones. The 
symmetric and anti-symmetric states differ from one another because the symmetric 
properties of their wave functions are different. To the best of our knowledge, in the 
papers published by far, except of our short communications [211 [22]. the consequences 
of this essential difference between the symmetric properties of the electron's wave 
functions were not treated properly. It is worth to be mentioned also, that the 
magnetotransport phenomena in such structures were not analyzed previously by taking 
into account the accurately calculated energies of Landau levels (LL). 

In this paper the parallel magneto-transport in DQW-structures composed of two 
identical QWs separated by thin barrier is studied. The experimental results are 
obtained for two different structures, the first one included QWs of quasi-rectangular 
shape, while the second one included the triangular QWs. Both of the structures are 
based on the InGaAs/InAlAs-system and for both of them we observed the beating 
effect in SdH-oscillations in perpendicular magnetic field arrangement (that is, without 
B\\ -component of the field). As it will be seen below, in order to interpret the beating 
effect which occurs in SdH-oscillations in the DQW with the matched density of states 
and in the absence of -component of magnetic field, it is necessary to introduce two 
different quasi-Fermi levels for the symmetric and anti-symmetric electron sub-systems. 

2. Description of the structures 

The DQW based on InGaAs/InAlAs/InP structures were produced by means of the low 
pressure metal organic vapour phase epitaxy (LP-MOVPE) on semi-insulating (100) 
InP: Fe substrates at the Institute of Electronic Materials Technology, Warsaw. The 
method of producing the structures with a single QW, which was reported earlier [23[ 124], 
was used also for producing DQW structure with rectangular, as well as triangular QWs. 
The structure of the first type (#2506) consists of two Ino.65Gao.35 As QW of about 20 
nm thickness each, and three Ino.52Alo.48As barriers. In each barrier there was the donor 
5-doping layer (see Fig. 1). 

The cross section of #2506-structure is depicted in Fig. 2a. The second DQW 
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Table 1. Parameters of the InGaAs/InAlAs DQW structure. 

Concentration Thickness Number of 5-doping layers 

Sample of (%) InAs in of the QW-profile (5-doping with donor concentration 

the channels channels layers (according to technologists) 



#2506 65 
#3181 65 



20 nm smooth interface 3 
5 nm sharp interface 2 



3.5 xlO^^cm-^ 
3.5 xl0^2cm-2 



In^Ali.^As 




In A 1, As 



In,Ga,.,As 



Ii\Al,_^As 



In Ga, As 




20 20 20 

Figure 1. The conduction band edge profile of the #2506 structure. 
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Figure 2. a) The cross-section of the #2506 structure used in the experiments; b) 
upper one: the piece of wafer with the indicated parts of it from where the #25061 
and #250611 was cut out; lower one: view of the sample; c) The cross-section of the 
#3181 structure. 



structure (#3181) whose cross-section is shown in Fig. 2c, also consists of two 
Ino.65Gao.35 As QW of the thickness of about 5 nm and three Ino.52Alo.48 As barriers. 
In this case however, the barrier between QWs was not doped. The corresponding 
conduction band shape is shown in Fig. 3. For this particular case the 2DEG in each 
QW could be considered as to be putted into a triangular potential well, since as it is 
shown in Fig. 3, the bottom of the conduction band in each QW makes a triangle with 
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Figure 3. The conduction band shape for the sample #3181. 



Both of the structures are symmetric with respect to the middle of central barriers 
this means that quantum wells are identical for both of them and with great probability 
their 2DEG densities are matched. The parameters of the structures of both types are 
listed in Table 1. 



3. Experimental results 

The magneto-transport measurements were performed by means of superconducting 
magnet, which gives the possibility to get the magnetic fields up to 11 T. The sample 
was mounted in the anti-cryostat which enables to change the temperature ranging from 
0.4 K to 300 K[2S]. The external gate voltage was not applied because this field destroys 
the symmetry of the DQWs profiles. The values of Hall {Uxy{B)) audi, transversal Uxx{B) 
voltage were recorded for the two opposite directions of magnetic and electric fields. 
Therefore, eight records (for magnetic field once going up and then going down) were 
made and averaged at each single measurement for the transverse magnetoresistance 
Rxx{B) and Hall resistance Rxy{B) at given temperature. These curves are shown in 
Figs. 4a-b and Figs. 5a-b. The pronounced SdH oscillations, which manifest themselves 
in the Rxx- magnetic field dependence, as well as Integer Quantum (IQHE) Hall effect 
(i?3,y(i?)-curves) are clearly seen. The SdH oscillations undergo also the beating effect 
which was observed in our experiments for both types of the structures. We observed 
four beating-effect nodes for the structure 7^25061 (see Fig.4a) and three nodes for 
the sample #3181 (Fig. 5a). It is important to mention that no beating-effect was 
observed in case of the structure with single QW based on the same heterostructures 
Ino.65Gao.35As/Ino.52Alo.48As[2l] . 

The several plateaux are observed in Rxy{B) curves for the magnetic fields up to 
6T for the sample #25061, as well as for #3181. It can be seen, that the plateaux 
observed in the magnetic fields greater than 6 T, are distorted. There could be two 
possible explanations of this distortion. The first one accounts for the presence of the 
additional conducting channels, for instance, due to the 6- layers. However, the effect 
of the additional conducting channels should be more pronounced at 4.2 K, rather than 
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Figure 4. a). The SdH curves for the structure #25061 at two temperatures: 0.6 
K and 4.2 K. b) The results of QHE measurements for the structure #25061 at the 
temperatures 0.6 K and 4.2 K. 
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Figure 5. a) The SdH curves for structure #3181 at 4.2 K. b) The resuhs of QHE 
measurements for structure #3181 at 4.2 K. 



at 0.6 K, while as it is easily seen, the minima of the distorted plateaux are deeper 
at 0.6 K than those at 4.2 K (see Fig. 4b). The second explanation of the plateau 
distortion refers to the QWs coupling. The additional maximum on the plateaux was 
already observed earlier (see Ref. 7). 

In Fig. 6 the results of another series of measurements of QHE and SdH oscillations 
are presented. These results were obtained also for the structure #2506, but the sample 
however was cut out from another piece of wafer (see Fig. 2b), that is why we denote it 
as #2506H. It is clearly seen that in this case the QHE-plateaux are flat and the Hall 
resistance Rxyis quantized at h/e'^ where the occupation factor is z/ = {No + Ni)sD. Here 
Nq is the uppermost Landau level occupied by electrons in the lowest subband (i = 0), 
A^^i means the same for the next subband (i = 1), s = 2 is the spin degeneracy, D = 2 is 
the symmetric degeneracy (symmetric and anti-symmetric states). At the same time it 
is necessary to note that together with regularly flat plateaux in QHE for the magnetic 
field ranging from 3 T up to 8 T, at the magnetic field 8.5 T the small distortion of 
the plateau is observed. We do not observe the beating effect in SdH oscillations for 
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the sample #250611 (see Fig. 6). In the next sections we will prove that the beating 
effect in SdH oscillations is caused by SAS-gap and the difference between two samples 
of structure ^^2506 is related to the difference in carrier density. 
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Figure 6. The SdH oscillations and QHE for the second sample #250611 of the 
structure #2056. 



The barrier width in the structure #2506 is equal to 20 nm and the initial, that is 
in zero-magnetic field, SAS-gap is so small, that can be completely screened by 2DEG 
(see Ref. 10 and Sect. 5 of this paper) in case of the sample #250611. In case of 
the sample #25061 the screening is insufficient to cause the collapse of SAS-gap. It is 
even more valid for the sample #3181, where the barrier width is 5 nm and SAS-gap is 
significantly larger. That means that in our case the initial SAS-gap for the two samples 
of the same #2506-structure can fluctuate because the electron density slightly changes 
along the radius of the wafer. Indeed, the electron densities determined by the slop of 
Rxy{B)-cuTve in a small magnetic field for the #250611 sample is greater then for the 
#25061. It is equal to 4.2 xlO^^ ^m-^ for the sample #250611 and 3.2 xlO^^ ^,^^-2 
#25061. As it is shown in Fig. 2c, the distance between two contacts for Vxx{B) is about 
0,8 mm, while the distance between two samples: #25061 and #250611 measured along 
the radius of wafer is 25 mm. As it follows from our experimental data, the difference 
in carrier concentrations for two samples #25061 and #250611 is about 30%. Since 
the distance between the tips of probes measured along the radius of wafer is about 0.8 
mm, the carrier concentration change within each of the samples is less than 1%. Hence, 
the homogeneity of carrier density is quite satisfactory within the samples used in our 
experiments. So, for the sample #25061 we observe at the same time the beating effect 
and distortion of the plateaux in QHE, while for the sample #250611 the beating effect 
is absent and the plateaux are fiat within the experimental error for the magnetic field 
< 8T. 
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4. Simulation of beating effect 

We carried out the Fourier analysis of the SdH oscillations for both structures, #25061 
and #250611, as well as for #3181. It is clearly seen that in case of sample #25061 
(Fig. 7a) there is predominating main harmonic which is splitted. The similar situation 
is observed for the sample #3181 (Fig. 7c) but in this case the magnitude of splitting 
is larger. That corresponds to beatting effect observed for these two samples. For the 
sample #250611 (Fig. 7b) there is one strong harmonic too but it is not so dominate 
as in previous two cases. It is possible to assume that observed number of harmonics 
suggests that no less than two subband are occupated by electrons in case of the sample 
#2506. 




Figure 7. The SdH-oscillations' Fourier transform plotted versus Bf — [A(1/_B) ^]: 
a) for the #25061 structure, b) for the #250611 structure and c) for the #3181 structure. 



From the other hand, as it was mentioned already, our earlier experiments with the 
structures containing a single quantum well[2l], did not show beating effect which could 
be attributed to spin splitting. Therefore we can conclude, that another kind of electron 
states splitting at zero magnetic field also exists. In order to reveal this most important 
effect, we have used the next method of the experimental SdH curve simulation. The 
starting point for further analysis is the next expression, which is well-known in the 
theory of SdH [26]: 



Ap _ ^ 5 / IP (3Tm' cosjlTfu) 
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exp{-l (^Tom /B) cos 2tt{1/PB -1/8- h). 



(1) 



This formula takes into account the unparabohcity of conduction band as well as 
spin splitting. Here Apis the deviation of p from background resistivity, po is zero-field 
resistivity, B is the transverse magnetic field, T is the temperature, f3 = 2TT'^kBmQ/he, 
P = he/Efm* stands for the SdH period and, is the Dingle temperature. 

The simulation procedure was based on adding only two Fourier harmonics of 
nearly equal frequency. According to that was mentioned above, the formula (1) can be 
simplified and reduced to two entries of the Fourier series only. As a result we have: 



and 



exp — — cos 



2 X {uji 



kB) X I 



exp — — cos 



TT 



2 X (cu2 + kB) X — 
B 



(2) 



(3) 



where: 71 ,72 are the coefficients which are equivalent to Rj3Tj:,m' in Eq. (1), cui, 
U2 are the cyclotron frequencies for two sub-systems of Landau levels. The kB-teim. 
is included into the argument of cosine in Eqs. (2), (3) in order to have an additional 
degree of freedom to get better agreement between experimental and calculated curves. 
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Figure 8. a) The upper curve represents the experimental data for SdH oscillations 
(#25061). The lower one represents the results of simulations performed according to 
the method described in Sect. 4. b) The same is for the structure #3181. 



It turns out, that in order to get better agreement between the experimental results 
and numerical calculations, one should reduce the oscillation frequencies when the 
magnetic field increases. These one could confirm the experimental fact that SAS- 
gap increases if magnetic field increases and if one suppose that zero magnetic SAS-gap 
causes the beating-effect. 

The results of the simulation are shown in Figs. 8a and 8b for the samples 7^25061 
and #3181, respectively. The fitting parameters used in the simulations, were chosen to 
be: 71 = 8.8, 72 = 8, Ui = 34.9, 002 = 40.27, k = 0.095, and 71 = 9, 71 = 8.8, ui = 30.8, 
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UJ2 = 23.8, k = 0.095 for these two samples, respectively. In this way we have obtained 
two oscillation components for each SdH-curve. These curves are shown in Figs. 9a, 9b 
and 10a, 10b. They should be interpret using of Landau levels energy calculates carried 
out for the symmetric and anti-symmetric states separately, by means of the best fit 
procedure. 

5. Calculation of Landau-level energy 

In our calculations we have used the model proposed in Ref. 10, where the symmetric 
DQWs separated by the middle barrier were considered. An external magnetic field 
was applied perpendicular to the planes of QWs. The electron's behavior was described 
by the Schrodinger equation which included the DQW potential, as well as the self 
consistent Hartree potentials. The in-plane movement in each of the subbands of the 
DQW was supposed to be Landau quantized. The total electron energy which includes 
the contributions of the in-plane as well as vertical electron motion, can be written as: 



where is the exchange energy (which involves screening) given by Eqs. (12)-(18) of 
Ref. [10], for DQW based on the GaAs well- material and Alo.sGao.yAs barrier- material 
with the well width of 14 nm and the barrier width of 3 nm. The authors of Ref. 10 
have found an approximately linear B-dependence of the tunneling gap for the magnetic 
fields < i? < 9T. One can conclude therefore, that the total splitting of energy states 
of two QW's At, consists of two parts: constant Asas, determined by the overlapping of 
the electron wave functions of two QW's, and changeable tunneling term V^^, involving 
screening factor: 



The last term, being strongly reduced when B increases up to the values for which 
the quantum limit occurs, can be written as: 



where Kqis the maximum screening factor in zero magnetic field depending mainly 
on the electron densities and —{kB)is the changeable screening factor decreasing with 
B up to the value when the first Landau level in two QWs becomes to be occupied. 
As it was observed experiment ally [TO] . the total energy gap Af between symmetric and 
anti-symmetric states in DQW is proportional to the magnetic field. 

We adapted this model for our particular case of symmetric Ino.65Gao.35As/Ino. 52 Alo.48As 
DQW structure (#25061) assuming rectangularity of QWs. The last one enables to cal- 
culate the Landau level energies according to the formula [27]: 




(4) 



A, = AsAS + Vni- 



(5) 




(i^o - kB) 



(6) 



(E'-E^)(E, + E' + E^) 




(7) 
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where 



/i 




■'^m* V 2/ 2 



c 



(8) 



E^ + Eg + |A 
(E, + A)(E^ + E, + |A) 



iEg + lA)iE^ + Eg + A) 

where E' = E ± (A^^^ + k' E±), Eg is the energy gap (between the top of valence band 
and the edge of conduction band) and is spin-orbital splitting, mo is the electron mass in 
a vacuum, m* is the effective electron mass corresponding to the edge of the conduction 
band, fis is the Bohr magneton, n = 0, 1,2... stands for Landau levels, i = 0, 1,2... is 
the number of sub-bands, E± is the Landau level energy for bulk semiconductor (we use 
here the W. Zawadzki's notation[27j) and E is the unknown n*'' Landau level's energy 
for the i*'^ sub-band in QW. The parameter A|^g is the energy gap determined by the 
overlapping of the wave functions, which is diminished by the screening at zero magnetic 
field as it follows from (5) and (6): 

At = {AsAs - Ko) + k'B = A'sAs + k'B. (10) 

Thus, the SAS-gap is proportional to the magnetic field B. The coefficient of 
proportionality k' can be considered as fitting parameter and is equal to 2k = 0.19, 
as it was determined in previous Section. We assume the triangular potential shape in 
case of asymmetric DQW (t^SISI) with one donor (5-layer in the barrier on one side 
of each QW. It means that Eq. (4) has to be modified for QW with the triangular 
potential (see Ref. 27): 



+ b)a'/^'/^ + {a-bY\n 



61/2 _ ^1/2 



(6-a)V2 



2m! 



1/2 

AeFhn{i + 3/4), (11) 



where: a = E' - E^; b = E*g + E' + E^; E' = E ± (A^^^ + k' E^); F is the electric 
field caused by the existence of the interface and which is determined by the potential 
U = eFz. The magnitude of E±_ is determined in the same way as previously, by means 
of Eqs. (8) and (9). The equations (7) and (11), together with (8) and (9) allow to 
calculate the Landau levels (LL's) energies in rectangular and triangular DQW with the 
spin-splitting and tunneling gaps variable in magnetic field taken into account. We used 
then the results of these calculations to account for the SdH oscillations and to find the 
parameters of 2DEG for both of the structures (the Fermi level (FL), carriers density 
etc.). 

The calculations of the LL's for the sample #25061 were performed using the band 
parameters, determined previously for the sample with the single quasi-rectangular 
QW[24]. These parameters are listed in Table 2, they are the same for both of the 
samples. As for the value of A|^^ = 1 meV for the #25061 structure and A|^^ = 5.5 
meV for the #3181 structure. 
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Table 2. Band-structure parameters used for the LL's energies calculatfon, where 
i?g-energy gap, TO*/mo-the ratio of electron mass and effective electron mass, A-spin- 
orbital splitting. 

X Eg,[eV] m*/mo A,[eV] 
0.65 0.723 0.0229 -9.90 0.355 



6. Discussion 

6.1. The role of SAS-gap. 

As was mention about since no beating-effect in SdH oscillations was observed in case 
of a structure with a single QW[23], we presume that this effect cannot be attributed 
to the spin-splitting of Landau levels. There are two alternative explanations of this 
effect. The first one is due to the existence of the SAS-gap at zero magnetic field. The 
second one is due to different carrier densities in two QWs (so called mismatching of 
carrier densities). It is well known that for different carrier densities in two QWs, the 
Fermi levels in both of them are also different if the tunneling effect is negligible. As a 
results of this difference, the beating effect could be observed. It is worthy to mention 
however, that according to Ref. 6 where the mismatching of electron densities were 
created in controllable way, the beating effect was not observed. We believe that in our 
case the beating effect is not the result of difference in carrier densities in two QWs, 
but is the result of the existence of SAS-gap in zero magnetic field. There are three 
evidences for that. The first one is the symmetry of the structures in questions with 
respect to the center of the middle barrier: both of the QWs were produced in a way 
to make them identical. The second one is that the width of the central barrier is small 
enough in order to equalize the carrier densities in both of the QWs of the structure 
#2506, as well as #3181. Note that in case of #3181-structure the barrier width is 
equal to 5 nm and in case of #2506-structure is 20 nm. If the beating effect were caused 
by the difference in electron densities, it would be more pronounced in case of #2506- 
structure then in #3181 and the period of beating (versus 1/B) would be smaller for 
the structure #2506, while in our experiment it was just the opposite. And at last the 
third evidence. In two samples of the same structure (#25061 and #250611) the electron 
densities are different, but the difference is comparatively small. If the beating effect 
were caused by the difference in electron densities in two QWs, it would be impossible to 
explain, why the beating effect disappears when the electron densities increase a bit in 
one of the samples (#250611). And the contrary, the disappearing of beating in #250611 
can be very naturally explained, if we attribute it to the existence of SAS-gap in zero 
magnetic field, because the small increase in electron density leads to the collapse of 
small SAS-gap due to screening. 

Therefore, the non-zero SAS-gap causes the SdH- oscillations at two different 
frequencies (versus 1/-B) and leads to beating effect. 
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In the next sub-section we will show that our experimental data related to beating 
effect can be sufficiently well interpreted, if one takes into account the accurately 
calculated LL energies and introduces two quasi-Fermi levels characterizing two electron 
sub-systems which belong to symmetric and anti-symmetric states in DQW-structure. 

6.2. Two quasi-Fermi levels. 

The calculations of LL's energies were performed according Eq. (7) for the sample 
#25061 and Eq. (11) for the sample #3181 using parameters shown in table 2. The 
results of our calculations for the sample #25061 are shown in Fig. 9a (for the symmetric 
states) and in Fig. 9b for the anti-symmetric ones. The calculated energies are plotted 
versus magnetic field for two subbands (i=l,2). Usually the Fermi level (FL) position are 
determined by LL's and FL crossings corresponding to the maxima of two components 
of SdH oscillations, acquired by simulation of beating effect. It is worth mentioning 
that, as it is seen in the Figs. 9a, 9b, the FL-positions are perfectly well fitted to the 
oscillation curves and are unique, since they are correspond to the regularly spaced 
crossings of LL's, equidistant with respect to the inverse magnetic field. 




Magnetic Field, T Magnetic Field, T 

Figure 9. These figures represent the best fit of the Fermi level energy corresponding 
to: a) symmetric part of the SdH oscillations b) anti-symmetric part of the SdH 
oscillations, in case of the structure #25061. 

The FL displacement up or down, destroys this regularity immediately, due to the 
superposition of the two space-quantized subbands. We should admit that we could not 
achieve the satisfactory agreement between our experimental data and the theoretical 
predictions about the SdH oscillation peaks positions, if we tried to use the single Fermi 
level common for both types of the states, symmetric and anti-symmetric ones. As 
it is clearly seen in Fig. 9a and 9b, the Fermi levels for two oscillation components 
are different, and the difference is equal to 12 meV. The similar is true for the sample 
#3181. In case of this sample only one subband is occupied, therefore to determine FL 
is simpler problem to solve. In spite of this, it was impossible to find a single Fermi level 
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Figure 10. These figures represent the best fit of the Fermi level energy corresponding 
to: a) symmetric part of the SdH oscillations b) anti-symmetric part of the SdH 
oscillations, in case of the structure #3181. 

which would be common for both sub-systems, symmetric and anti-symmetric. It turns 
out that the difference between two Fermi levels for the symmetric and anti-symmetric 
states in this case is equal to 32 meV (see Fig. 10a and 10b, respectively). 

It is interesting to note, that in case of the sample #25061, the value of total splitting 
At is approximately equal to 15 meV for n=0, i=0, in a magnetic field of about B=10 T, 
and should increase if the field increases up to about 30 T, where the predicted quantum 
limit should occur. The obtained values of the FL energies give us the possibility to 
estimate the density of carriers for both structures. These values are 3.8x10^^ cm~^ for 
the sample #25061 (good agreement with date obtained by the slop of Rxy{B)-c\XT've in 
small magnetic field) and 2.5x10^^ cm~^ for the sample #3181. In the first case this 
value is a bit greater than the technologists set it up to be, and in the second case is a bit 
smaller. So, our experimental data force us to introduce the concept of two quasi-Fermi 
levels, by means of which two electron sub-systems can be characterized. These two 
sub-systems are the electrons belonging to two different types of states, symmetric and 
anti-symmetric ones. One can encounter the concept of quasi-Fermi levels in different 
contexts for instance, in semiconductor physics. The quasi-Fermi levels for electrons and 
holes very often are used; quasi-Fermi levels for the electron states with and k^{kaie 
for the electron wave vectors along say, +x and —x axes) are introduced to explain the 
peculiarities of IQHE|28[ 129]. and so on. Generally speaking this concept simply means, 
that one can treat two sub-systems as weakly interacting and having their own rate of 
establishing the equilibrium state. Therefore, our next step is to argue that indeed, we 
could interpret our data, introducing the two quasi-Fermi levels which characterize two 
separate electron sub-systems. 
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6.3. Electron- electron interaction for the two sub-systems: symmetric and 
anti- symmetric states. 

Suppose now that the 2DEG in the DQW structure under consideration is off the 
thermodynamic equihbrium state and the deviations from the equihbrium are small 
enough. Obviously, these deviations are caused by the current which flows through the 
structure. Then the electron distribution functions, one of which corresponds to the 
symmetric state f^""^ and the other to anti-symmetric one f^''\ can be represented as 
follows: 

where 

Here the superscript a stand for the symmetric and b for the anti-symmetric states (see 
App. A where the wave function for this states are considered). Denote the electrons 
belonging to sub-system a by means of ipa, and belonging to sub-system b as ipb and 
define the rate of establishing the equilibrium state in each sub-system as: 
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are the probabilities for the electrons belonging to the sub- 



system a or 6 to be scattered from the state kto the state k'and vice versa (see Fig. 11). 
These probabilities are determined by the next formula: 
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E^^'''^ are their energies 



are the initial and final states, E^"'''^^ 
and V is the operator, corresponding to the perturbation responsible for the transition. 
It is well-known that mainly the electron-electron (e— e) collisions lead to establishing the 
equilibrium states within the electron gas in a semiconductor at low temperatures [30] . 
It means that whatever the initial distribution function is, the final distribution can be 
described by displaced Fermi (or displaced Maxwellian) distribution function. However, 
for establishing of such distribution, an important condition has to be satisfied: e — e 
collisions should be more frequent than the electron-phonon (e — ph) collisions. In other 
words. Tee << T^-ph, whcrc Tee and Te-ph are the e — e and e — ph scattering times, 
respectively. Since the temperatures at which our experiments were carried out, were 
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very low, we can suppose the last condition to be fulfilled. Now in order to calculate 
Tee = T^"''''^k foT the two sub-systems, we have to calculate W^°^l, the corresponding 
transition probabilities. Therefore, we should analyze the probabilities for two electrons 
to transit from the states Uifci,/2^2) to the states [l^k^.l^k^) where Iil[,l2,l2 
stand for the quantum numbers, other then k. 




k 



Figure 11. Diagram, representing the electron-electron scattering in case of two sub- 
bunds, symmetric and anti-symmetric. 



Notice that the subscripts 1,2 stand for the electrons 1 and 2 in a pair, while 
the superscripts (a,b) are to distinguish symmetric and anti-symmetric states of the 
electrons in DQW structure. Then 

X 5 (e^n e^"^-") - 6>'^) - e^")) . (16) 

Here we introduced the following notations liki = 1, l[k[ = l' etc; e^"'^'', €2"'^'', ei"''^\ 
€2 "''''^ , are the electron energies in the initial and final states, respectively. It is necessary 
to calculate the matrix element M^'^2^^(^i, 2/) In Eg- 16 to show the different rates of 
establishing of the quasi-equilibrium states in two sub-systems. The calculation of two 

elements ^^^"2)^(1' 2') ^('1 2)^(1' 2') caring out in App. B. 

By inspecting Eqs. (B.3) and (B.4), even without doing out the evaluation of these 
integrals, one can easily see that M^^i 2)^(1' 2') 7^ ^{i2)-*{i' 2') ^^"^ hence, the same is true 
for Tg^e'*-*: r]g) 7^ r^g^. That is, the rates of establishing of the quasi-equilibrium state 
under electron-electron scattering in two sub-systems, corresponding to the symmetric 
and anti-symmetric subbands, are different. 

The same arguments are also valid, if we consider the electron-phonon interaction. 
Moreover in case of electron-phonon interaction the formulae analogous to (B.3) and 
(B.4) are even simpler, because in this case the initial as well the final electron states 
are one-particle states, but not the two-particles ones, as in case of electron-electron 
scattering. Anyway, one can state that both rjg'''^ and Tg"p)i in two sub-systems 
corresponding to the symmetric and anti-symmetric states, are different: Tg^g-* 7^ rjg^; 

e—ph I e—ph " 
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The only thing which also has to be proven, is that these two sub-systems are weakly 
interacting. To this end let us consider the transition matrix element ^(^i 2)^(1' 2') where 
(1,2) belong to initial symmetric state ip'} and (1', 2') - to final anti-symmetric state ■?/'/• 
Then, under the integral sign in (B.3) one should have the product of two functions 
ipa{,zi), Lpi,{zi)^ the one is even and another one is odd. Notice also, that U'^^is an even 
function of their variables. Hence, M^^2)^{v 2') should be very nearly to zero due to 
parity selection rule, and it means that the electron-electron scattering does not mix the 
states of the different parities. Of course, electron-phonon interaction can and usually 
do them mixed. However, the temperatures at which our experiments were carried out 
were so low and the electron densities in QW so high, that we the main scattering 
mechanism is the electron-electron scattering, and that two sub-systems belonging to 
symmetric and anti-symmetric states are weakly interacting in that sense. 



7. Summary 

We have studied the parallel magneto-transport in DQW-structures of two different 
potential shapes (quasi-rectangular and quasi-triangular) and have found an important 
contribution of the symmetric properties of the charge carriers eigenstates to the 
magneto-transport phenomena. The beating effect occurred in SdH oscillations was 
observed for both types of DQW. To explain this effect, we developed a special scheme 
for the Landau levels energy calculations based essentially on the model proposed in 
Ref. [10] and carried out the necessary simulations of beating effect. In order to obtain 
the agreement between our experimental data and the theoretical predictions concerning 
SdH-oscillations, we introduced two different quasi-Fermi levels, which characterize sub- 
systems of symmetric and anti-symmetric states in DQWs. The existence of two different 
quasi Fermi-levels simply means, that in DQWs consisted of two identical InGaAs 
QWs separated by InAlAs-barriers no wider than 20 nm, the total electron system can 
be considered as divided into two sub-systems characterized by symmetric and anti- 
symmetric wave functions respectively. One can treat then these two sub-systems as 
weakly interacting in the sense that electron-electron scattering does not mix them and 
that they are described by their own distribution functions characterized by their own 
quasi-Fermi levels. 



Appendix A 



The quantum magnetotransport in DQW can be treated in the framework of the effective 
mass approximation [29j. The main equation of the effective mass model is of the form: 

, 2 



(ihv 



eA 



2m* 



+ Uiz) 



il){r) = Eip{f) 



(A.l) 



where Eg = Ec + ei, E^. stands for the bottom of the conduction band, Si is the 
energy of the first subband due to space quantization along 2;— axis and U{z) stands for 
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the potential shape along z-axis, as in case of Fig. 1 or Fig. 3. In order to take into 
account the magnetic field arrangement e^is the unit vector of ^— axis, and BJ^j, 

where j is the current density) choose the next vector potential gauge: = By, Ay = 0, 
= 0.) Look again at Fig. 1 and Fig. 3; what is peculiar about these pictures is that 
regardless of the shape of QWs, rectangular or triangular, the structures are symmetric 
with respect to the centre of the barrier which is in the middle of the structure. The 
Hamiltonian, which describes the electron moving in the potential like that of Fig. 1 
and Fig. 3, is invariant under inversion z — > —z and hence, its eigenfunctions are 
double degenerate. However, due to the tunneling effect this degeneracy is removed and 
the lowest energy state in each quantum well becomes splitted into two, bounding and 
anti-bounding states represented by even and odd functions, respectively. We termed 
these states also as symmetric and anti-symmetric ones. Then the eigenfunctions of the 
Hamiltonian (A.l) are of the form: 

V'a(6)(r) = exp{ik^x)Xm{y)Va{b)iz)- (A.2) 

In the last expression k^is the wave vector along axis, Xm{y) is the wave function 
corresponding to the m*^ Landau level, is the structure length along x— axis, ipa{b){z) 
and stands for the anti-symmetric (symmetric) state corresponding to the symmetric 
structure of the Fig. 1 or Fig. 3. 



Appendix B 

Two electrons interacting in the e—e scattering, are indistinguishable and hence, the two 
electron states (1,2) and (l',2'), before and after scattering event, should be described 
by the Hartree-Fock determinants as follows: 

V'l (rl) ^2 (r1) = [V'l (fl) (r-^2) - V'l (r,) (fl)] . (B.l) 
In order to calculate the matrix elements entering the expressions for W^'^'^'l i^, in 

k^k' [k'^kj 

(14) using the last expression and (A.2), one should substitute (fi{fi) <f2{r2), — 1, 2) 
by the next functions: 



V'jfe (n) = -jr= i^f^jx^i) Xmj ivi) Vb (zi) = 



1 



r ^ki i^i) Xmj iVi) Vb {Zi) , 



i^ja (rl) 



1 



: exp {ikj^Xi) Xmj iVi) fa {Zi 



--¥>ki i^i) Xmj iVi) fa (Zi) , 



(B.2) 



where ipa{zi) and ipb{zi) correspond to symmetric and anti-symmetric states, 
respectively. Here in the last expressions kj = kj^- is the wave vector corresponding to 
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the electron movement along ,T~axis perpendicular to axis. Then the corresponding 
matrix elements -^(i 2)^(1' 2') form: 



''^(l,2)^(l',2') 

^/Tr I I ^^2V^fc'l(^l)^™'l^^l)'^^(^l)^fe'2(^2)X^2(^2)</:'6(^2) X 

U^^^k2{Xl)Xm2{yi)M^l)'Pkli^2)Xml{y2)M^2), (B.3) 

(1,2)^(1', 2') 

7^ / "^^^ / ^^2<^I-'l(^l)x;^'l(?/l)<^:(^l)¥'J2(''2)X^2(?/2)</':(^2) X 
C/^Vfel(a;i)Xml(yi)V'a(2l)V^fc2(^2)Xm2(y2)V^a(22) " 

dl^i y (il^2<^J/^l)x;'i(Z/l)<^:(^l)¥'j2(^2)x;'2(Z/2)¥':(^2) X 
U^^Vk2i^l)Xm2{yi)^aizi)(Pkiix2)Xml{y2)^aiz2), (B.4) 

where = dxidyidzi and (iV2 = dx2dy2dz2. The first term in each of these formulae 
corresponds to the direct interaction, while the second one to the exchange interaction, 
and 



jjee ^ jjee (^^^^ ^ 

■ (B.5) 

/cj y {xi - X2) + (1/1 - 1/2) + (^1 - Z2) 

Here £20 eg (^) is the dielectric constant of the two dimensional electron gas which 
incorporates all the screening effects in 2DEG. The integration in the formulae (B.3) and 
(B.4) is carrying out in such a way that at least one of the entries in denominator in (B.5) 
((,Ti — X2)'^, {yi — 1/2)^ or (-21 — Z2Y) is not equal zero, guarantying the convergence of the 
integrals. This is the consequence of the fact that Lpji,{fi), Lpja{fi) arc the wave functions 
corresponding to some definite energy and this energy determines the minimum distance 
at which two electrons can approach each other. 



References 

[1] M. C. Bronsager, K. Flensberg, B. Yu-Kuang Hu, Antti-Pekka Jauho, Phys. Rev. Lett. 77, 1366 
(1996) 

[2] H. Rubcl, A. Fischer, W. Dietschc, K. von Klitzing, K. Ebcrl, Phys. Rev. Lett. 78, 1763 (1996) 
[3] X. G. Feng, S. Zclakicwicz, H. Noh, T. J. Ragucci, T. J. Gramila, Phys. Rev. Lett. 81, 3219 (1998) 
[4] J. G. S. Lok, S. Kraus, M. Pohlt, W. Dietsche, K. von Klitzing, W. Wegscheider, M. Bichler, Phys. 
Rev. B. 63, 041305(R) (2001) 



Magnetospectroscopy of symmetric and anti-symmetric states in double quantum weUs20 

[5] M. Kellog, I. B. Spielman, J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 88, 
126804 (2002) 

[6] H. S. Fresser, H. Frey, F. E. Prins, D. A. Haram, D. P. Kern, J. Bottcher and H. Kunzel, Semicon. 
Sci. Technol. 15, 242 (2000) 

[7] G. S. Boebinger, H. W. Jiang, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 64, 1793 (1990) 
[8] K. Ensslin, A. Wixforth, M. Sundaram, P. F. Hopkins, J. H. English, A. C. Gossard, at. all., Phys. 
Rev. B. 47, 1366 (1992). 

[9] K. M. Brown, N. Turner, J. T. NichoUs, E. H. Linfield, M. Pepper, D. A. Ritchie, G. A. C. Jones, 
Phys. Rev. B 50, 15465 (1994) 
[10] Danhong Huang, M. O. Manasreh, Phys. Rev. B 54, 2044 (1996) 
[11] J. P. Eisenstein, A. H. MacDonald, Nature, 432, p.691, 9 December 2004 

[12] R. M Lewis. P. D. Ye, L. W. Engel, D. C. Tsui, L. N. Pfeiffer, K. W. West, Phys. Rev. Lett, 89 
136804-1, 2002 

[13] R. M Lewis, Yong P. Chen, L. W. Engel, D. C. Tsui, L. N. Pfeiffer, K. W. West, Phys. Rev. B., 
71 081301(R), 2005 

[14] G. Gervais, H. L. Stormer, D. C. Tsui, L. W. Engel, P. L. Kuhns, W. G. Moulton, A. P. Reyes, L. 

N. Pfeiffer, K. W. Baldwin, K. W. West, Phys. Rev. B., 72, 041310(R), 2005 
[15] L. A. Tracy, J. P. Eisenstein, L. N. Pfeiffer, K. W. West, Phys. Cond. Mat. 0511321, 2006 
[16] G. S. Boebinger, A. Passner, L. N. Pfeiffer and K. W. West, Phys. Rev. B. 43, 12673 (1991) 
[17] T. Jungwarth, T. S. Lay, L. Smrcka, M. Shayegan, Phys. Rev. B, 56, 1029 (1997) 
[18] Y. W. Suen, J. Jo, M. B. Santos, L. W. Engel, S. W. Hwang and M. Shayegan, Phys. Rev.B 44, 

5947 (1991) 

[19] J. A. Simmons, S. K. Lyo, N. E. Harff and J. F. Klem, Phys. Rev. Lett. 73, 2256 (1994) 
[20] T. Ihn, H. Carmona, P. C. Main, L. Eaves and M. Henini, Phys. Rev.B 54, R2315 (1996) 
[21] E. M. Sheregii, D. Ploch, M. Marchewka, G. Tomaka, A. Kolek, A. Stadler, K. Mleczko, D. Zak, 
W. Strupinski, A. Jasik, R. Jakiela, Inst. Physics Conf. Ser., No 187, Paper presented at 12th 
Int. Conf. Narrow Semiconductors, Toulouse, France, 3-7, July 2005, Taylor and Francis, New 
York, 2006, pp. 585- 590 

[22] M. Marchewka, E. M. Sheregii, D. Ploch, I. Tralle, G. Tomaka, M. Furdak, XVIII Int. Conf. High 

magnetic fields in the Semiconductor Physics, Wurtzbourg 2006, (accepted for publication) 
[23] E. M. Sheregii, D. Ploch, M. Marchewka, G. Tomaka, A. Kolek, A. Stadler, K. Mleczko, W. 

Strupinski, A. Jasik, R. Jakiela, Parallel magneto-transport in multiple quantum well structures, 

Fizika Nizkikh Temperatur (English version: Low Temperature Physics edited by the American 

Institute of Physics) 30, no.lL 2004 pp. 1146-1156 
[24] G. Tomaka, E. M. Sheregii, T. Kakol, W. Strupinski, R. Jakiela, A. Kolek, A. Stadler, K. Mleczko, 

Crys. Res. Technol. 38, No. 3-5, 407-415 (2003) 
[25] The measurements at the temperature of about of 4.2 K were performed at Rzeszow University of 

Technology. 

[26] D. G. Seller and A. E. Stephens, /in: Landau level Spectroscopy, vol. 27.2 of series Mod. Problems 
in Cond. Matt. Science, Eds: G.Landwehr and E. I. Rashba, cds of series: V. M. Agranovich, 
A. A. Maradudin, Nort-HoUand, Amsterdam, pp. 1035- 1137 (1991) 

[27] W. Zawadzki, J. Phys. C 16, 229 (1983) 

[28] J. H. Davies, The Physics of Low-dimensional Semiconductors, Cambr. Univ. Press, Cambr. (1997) 

[29] S. Datta, Electronic Transport in Mesoscopic System, Cambr. Univ. Press, Cambr. (1999) 
[30] B. Gantmacher, Y. Levinson, Carrier scattering in metals and semiconductor, North-Holland, 
NY-Tokyo (1987). 



